LIMITING SET OF SECOND ORDER SPECTRA 



LYONELL BOULTON 

Abstract. Let M be a self-adjoint operator acting on a Hilbert 
space H. A complex number z is in the second order spectrum of 
M relative to a finite dimensional subspace C C DomM^, iff the 
truncation to C of (M — z)^ is not invertible. This definition was 
first introduced in ^ and according to the results of [J] and \V^\ . 
these set provide a method for estimating eigenvalues free from the 
problems of spectral pollution. In this paper we investigate various 
aspects related to the issue of approximation using second order 
spectra. Our main result shows that under fairly mild hypothesis 
on M, the uniform limit of these set, as C increases towards 71, con- 
tains the isolated eigenvalues of M of finite multiplicity. Therefore, 
unlike the majority of the standard methods, second order spectra 
combine non-pollution and approximation at a very high level of 
generality. 



1. Introduction 

Let M be a self-adjoint operator acting on an infinite dimensional 
Hilbert space Ti and let A be an isolated eigenvalue of M of finite 
multiplicity. Suppose we intend to estimate A numerically. A natural 
method, which is already over one hundred years old, is to truncate M 
by taking an orthonormal basis {(j)j}JLi of Ti contained in the (oper- 
ator or form) domain of M and compute the eigenvalues of the n x n 
matrix M„ = {{Mcpj, (l)k))]^k=i large n, expecting that some of these 
eigenvalues would be close to A. Unfortunately, an ill chosen {(pj}'^^ 
will result in either of the following difficulties: 

- Lack of approximation: no eigenvalue of M„ is close to A. This 
happens when M is unbounded. 
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- Spurious eigenvalues: or spectral pollution, M„ has many eigen- 
values unrelated to the spectrum of M in a large neighbourhood 
of A. This occurs when A is in a gap, i.e. in between two points, 
of the essential spectrum of M. 

It is well known that if M is an operator bounded below and 
A < min[SpeCess then the Rayleigh-Ritz theorem provides satis- 
factory procedures to deal with these two issues. For instance, suppose 
that A is the first eigenvalue of M and it is non-degenerate. If 
is such that 



then the first eigenvalue of M„ converges to A (i.e. we achieve approx- 
imation) and the second eigenvalue of M„, counting multiplicity, can 
not be smaller than the second eigenvalue of M (i.e. no chance of spu- 
rious eigenvalues). Condition (H)) is useful in applications, because it 
can sometimes be verified on abstract grounds without having much 
information about the eigenfunctions ip, cf. theorem XIII. 4]. 

In contrast, when A is in between two points of the essential spec- 
trum, it is well known that pollution can arise. In fact, for each fi G 
]R\Spec M lying between two points of the essential spectrum, there ex- 
ists an orthonormal basis and a subsequence n(l) < n(2) < . . ., 
such that /i G Spec Mn{k) for all G N (see j7]). There is a large amount 
of literature devoted to study this problem in applications, we refer to 
[S] and [7] , and the references therein, for a more complete account on 
the matter. 

In order to find eigenvalues in gaps of the essential spectrum, one can 
appeal to the following geometrically motivated idea (see [S] and 0). 
Suppose a <( <b, such that a, 6 G Spec ess-^ but (a, 6)nSpec ess^ = 0- 
Then (M — Q"^ is bounded below by zero and the eigenvalues be- 
low SpeCess(^ — C)^ are the (A — C)^ > 0, where A is an eigen- 
value of M inside (a, b). Hence by applying Rayleigh-Ritz principle to 
(M — (Y, we might be able to estimate (up to a square root ambigu- 
ity) the eigenvalue A. If, for instance, [M — satisfies (^, from the 
first eigenvalue of its matrix truncations we may estimate the closest 
eigenvalue of M to C in (a, b). 

Instead of looking for the eigenvalues close to zero for a truncation 
of (M — Q"^, we can also consider finding z such that the truncation 
of (M — is singular. Denote by £„ the Span{0j}"^]^ and by Pn 
the orthogonal projection onto £„. Assume that C DomM^. We 
say that 2; G C is an element of the second order spectrum of M 



(1) 
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relative to £„, denoted below by Spec2 (M|£„), if and only if 
det [P„(M — = 0. This definition was first introduced in j3] 

by Davies. Levitin and Shargorodsky have recently proposed, see j7], 
using second order spectra to approximate isolated eigenvalues inside 
spectral gaps, on the grounds that they are free from pollution. To be 
precise, if z e Speca (M|£„), then 

(2) [Rez- \lmz\,Rez+ \lmz\] (1 SpecM 7^ 

(cf. |7j and |13j). In other words, a point in Spec2 (M|£„) is close 
to the real axis only when it is also close to the spectrum. Therefore 
numerical procedures based on second order spectrum, never lead to 
spurious spectrum. 

Although the bound (j2I) does not guarantee that any point at all 
in the spectrum is estimated for large n, some numerical experiments 
performed in [7j and the results of |2], indicate that approximation 
occurs in test models. 

In this paper we discuss rigorously various aspects of the question of 
approximation in the method considered by Levitin and Shargorodsky. 
For this, we study the uniform limiting set 

A := u— lim Spec2 (M|£„,) = {z E C : z„ — > z, 2;„ G Spec2 (M|£„)}. 

n— »oo 

The interest in A is clear, the points in Spec M which are estimated by 
Spec2 (M|£„) for n large, are those in AflM. The following assumption 
on A and P„ will be crucial in our subsequent analysis, 

, . ifMip = then ||P„MP„V^ - A^|| ^ 

^ ' and WPnM'^Pnip - X'^^W ^0, as n^oo. 

Our principal task will be to show the following 

Theorem 1. If X is an isolated eigenvalue of finite multiplicity of M 
and (H) holds, then A G A. 

Obviously any bounded operator satisfies (H), and notice the analogy 
between (H) and (HJ in the general case. We will discuss in section 5 
how to verify this condition for unbounded operators without having 
much information about ip. 

We should emphasize the relevance of this result. In many appli- 
cations, the essential spectrum can be found analytically, whereas the 
isolated eigenvalues of finite multiplicity are the ones that should be 
estimated numerically. Theorem ^ shows that second order relative 
spectra provide a method which will give convergence to these eigen- 
values. Furthermore, in conjunction with (0), it is ensured that the 
method will not pollute, no matter the location of these eigenvalues. 
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We are not aware of any other technique that combines approximation 
and non-pollution at such level of generality. 

In section 2 we keep close to the ideas of 0. We characterize 
Spec2 {M\Cn) and the set A in terms of resolvent norms. This charac- 
terization will be important in section 3, where we prove theorem ^ 
In section 4 we discuss in detail a simple model bounded operator, for 
which A can be found explicitly. In order to test theorem^ in concrete 
situations, we also include some experimental numerical outputs for 
rank one perturbations of this model. Finally section 5 is devoted to 
discuss the verification of condition (H) for concrete unbounded situa- 
tions. 

2. Second order spectrum and norm of inverses 

The sequence of subspaces £„ does not necessarily need to be nested 
in order to accomplish approximation. Below P„ denotes a family 
of orthogonal projections and := RanP„ C DomM^, such that 
dim£„ = n and -P„x — > x for all x E Ti. For the linear operator iJ, 
we shall often denote by if„ the truncation PnH\Cn- By we shall 
always mean the inverse of Hn within £„. If !!„ is not invertible we 
allow the abuse of notation = oo- 

From the definition it is easy to deduce that Spec2 (M|£„) is a set of 
at most 2n different complex numbers. In general these points do not 
intersect the real line, unless contains an eigenfunction of M. Since 



Spec2 (M|£„) is symmetric with respect to M. The limiting set A might 
or might not contain non-real points. The reference [2] is largely de- 
voted to showing that A C Spec M, when M is a discrete Schrodinger 
operator whose potential decays fast at infinity. On the other hand, in 
the example we discuss in section 4, A is equal to the unit circle. 

In order to describe properties of A, it is useful characterize the points 
in Spec2 (M|£„) as the zeros of a certain Lipschitz continuous function 
an{z). Davies in and Davies and Plum in study a procedure 
for computing eigenvalues of bounded self-adjoint operators, based on 
evaluating an{z) at 2; G M. Let 



liz ^ Spec2 (M|£„), then a^iz)-^ = ||[P„(M - zf\Cn]'^\\. Moreover, 
an{z) = if and only if ^ G Spec2 (M|£„). Both assertions are direct 
consequence of the definition. 



det Pn{M - zy\Cn = det Pn{M -zY\C, 




The following result shows that, in order to find the zeros of (T„, 
we can start from an initial guess and move in steps towards a local 
minimum. 

Lemma 2. The function (T„(z) is a non-negative Lipschitz continuous 
function of the variable z. Its only local minima are the points where 
it vanishes. 

Proof. We follow closely [H theorem 21]. If |w — 2;| < e, then 

\\P4M - z)^\CnV\\ 
OrXz) = mf 

^ . ^ \\Pn{M - Wy\CnV\\ + \W-Z\ ||P„(2M -Z- w)\CnV\\ 

< mf 



(3) - M 

\\{2Mn - z-w)v\ 

< 0-n(.w) + \w — z\c, 



< 0-n{w) + \w — z\ sup 



where c > can be chosen independent of e. This ensures the Lipschitz 
continuity of a^. 

On the other hand, given zq E C \ Spec2 (M|£„), choose u,v E Cn 
vectors of norm 1 such that 

an{zo)-' = \ {[PniM - Zof\Cnr'u,v)\. 

The function ([P„(M— t>) is analytic in 2; G C\Spec2 (M|£„). 
We claim that it is not constant. If such a claim is true, then there 
always exists z in any neighbourhood of Zq such that 

aniz)-' > \ {[PniM - zf\Cn]-'u,v)\ 

> \ {[P^iM-zof\Cr,]~\v)\ 

= CTniZoy^, 

ensuring the second assertion of the lemma. In order to prove our 
claim, first notice that 

\{[PniM - zY\Cn]-'u,v)\ < \\[P^{M - zf\Cn]-'\\ 

= \z\-^\\[z-'^PnM'^\Cn - 2z-^PnM\Cn + 1]'% 

Then, if \z\ > 3max{||P„M2|/:„||, ||P„M|/:„||, 1}, 

\\z~^PnM^\Crr - 2z-'^PnM\CJ < 7/9 < 1 



and hence 



\{[Pr,{M-Zy\Cnr'u,v)\ < (9/2)|z|-2^0 
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as \z\ oo. Thus necessarily the abovementioned function is not 
constant. 

In conjunction with A, it is natural to consider the set of asymptotic 
zeros of the family {o"„}, 

S = S(M) := {2 G C : an{z) ^ 0, n ^ 00}. 

If M is bounded, this set contains the hmiting set A. Indeed, suppose 
that z ^ S. Then there is a subsequence n(j) and c > 0, such that 
P„,(j)(M — zY\Cn{j) are invertible for all j G N and 

c > a~l^p) = \\[PnU){M - zf\C^^,)]-% 

By virtue of (jHI), for 1^; — w| < e, 

(^n{j){z) < (Tn{j)iw) + \w - z\Ci 

where Ci > can be chosen independent of e and j (here we use that 
M is bounded in order to ensure that Ci does not depend on j). By 
letting e be small enough, we can find C2 > such that cr„(j)(w) > C2 
whenever |w7 — 2;| < e. Thus 

00 

U Spec2 (M|£„(,)) n{\z-w\<e} = 
i=i 

so therefore z ^ A. 

Furthermore, whenever M is bounded, S fl M = SpecM. Indeed, 
if A G SpecM, for each A; > there is ipk ^ 'H, Wi'kW = 1, such that 
\\{M - X)'^iljk\\ < l/k. Then 

n^oo \\Pn^k\\ 

and so crn(A) -^0 as n ^ cxo. Conversely, notice that 

(4) {Pn{M-z)^\CnV,v) = {{M-Zfv,v) 

for all V G £„, ||f || = 1. Then, if A G M but A ^ SpecM, 

NumP„(M- A)2|£„ c Num(M- A)^ C [6,00) 

for some constant 6 > independent of n. Hence 

a„(A) > dist [0, Num P„(M - \f\Cn] > b, 

so that A ^ E n M. This observation is crucial in the method of Davies 
and Plum mentioned earlier. Notice that the direction E fl M C Spec M 
does not require M to be bounded. Here and below "Num" denotes 
numerical range. 
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It would be of obvious interest to find hypothesis which guarantee 
A = S. We show that the equahty holds at least in particular cases. 
Below and elsewhere we denote by T the unit circle. 

Proposition 3. If Spec M consists of two points a <b, then 

^ = AC{zeC:\z-{a + b)/2\ = {b - a)/2}. 

Proof. Since 

Spec2 (aM + (3\Cn) = aSpec^ (M|/:„) + (3 

and S(q;M + P) = aT,{M) + /? for scalars a and /?, it is enough to show 
the desired property for the case SpecM = {±1}. By virtue of the 
hypothesis, it is clear that z = can neither be in S nor in any of the 
Spec2 (M|£„), so below we assume z 0. 

The latter assumption on M yields = Id. Then, since 

(5) P„(M - 1£„ = 2z[iz-^ + z)/2 - P„M|£„] , 

z G Spec2 {M\Cn) if and only if 

(z^^ + z)/2 e SpecM„ C NumM = [-1,1]. 

Hence all the second order relative spectra of M are contained in T. 
This shows the second relation. 

We already saw that A C S in general. In order to show the reverse 
inclusion, assume that z ^ A. Then there is a subsequence n{j) and 
S > 0, such that 

oo 

IJ Spec2 {M\CriU)) ^ {\^ -w\ <6} = 0. 
i=i 

By virtue of 

o-nij){z) = 2\z\ dist [(z^^ + z) / 2, Spec Mn{j)]. 
But if {w"^ + w)/2 is in the spectrum of Mn{j)^ 

\{z'^ -z)/2- {w-^ - w)/2\ = \z\^^\w - z\ \w^^ - z\ 

> \z\-'6\ 

Hence an{j){z) > c > for all j G N and so 2; ^ S. This shows 
proposition El 

In applications it is of interest to understand how the second or- 
der spectrum and its limiting set change under compact perturbations. 
The following lemma provides some lights in this respect. Let if be a 
bounded not necessarily self-adjoint operator and K he compact op- 
erator. It is well known that is uniformly bounded for all n 
large and H + K is invertible, if and only if \\{H + -ft')^^!] is uniformly 
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bounded for all n large and H is invertible (see, for instance ^ theorem 
2.16]). This allows us to prove that S does not change substantially 
under compact perturbation. 

Lemma 4. // M is bounded and K = K* is a compact operator, then 

S(M) U Spec disc(M + K) = S(M + K) U Spec discM. 
In particular S(M) and S(M + K) coincide outside the real axis. 
Proof. The point z is not in the set at the left hand side if and only if, 

{M + K- = {M- zf + [(M - z)K + K{M - z) + K^] 
is invertible and 

||[P„(,)(M-z)2|£„(,)]-i 

is uniformly bounded for some suitable subsequence n{j). According 
to the above observation, the latter is equivalent to having (M — z)"^ 
invertible and 

||[P„(,)(M + ir-;.)2|£„(,)]-ll| 

uniformly bounded, so therefore it is equivalent to z not being in the 
set at the right hand side. 

3. Proof of theorem 1 
We start by noticing that condition (H) ensures that 
(6) o"„(A) ^0 as n oo. 

Indeed, take v = Pni^, then 

|P„(M-A)2p„Vl| 



IIP.^II 



as n — >• oo. 
Decompose 

1-1 = 8® 8^ 

8 = Span{^ e DomM : Mip = \ip} C DomM^ 

where dim£^ < cxo, and both 8 and 8-^ are invariant under M. Accord- 
ing to the standard notion, 8''~ invariant under M means that Mx G 8'^ 
for all X G 8-^ fl DomM. The restriction M\8 corresponds to multipli- 
cation by A. Denote by Pg the orthogonal projection onto 8. Put 

M = M(/ - Pe) + (A - ix)Pe + /iP^ =: A + /iP^ 
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where 7^ /i ^ SpecM and DomA = DomM. Then A = A*, M is a 
finite rank perturbation of A and 

Spec A = {X- n}U Spec M \ {A}. 

For z eC, let A{z) := {A - zf with DomA{z) = DomM^. Then 
A(z) is a holomorphic family of type A for 2; G C. If K{z) = /i(2A — 
2z — fj,)P£, then 

(M - z^x = A{z)x + K{z)x, X G Dom M^. 

Put An{z) := PnA{z)\Cn and = PnK{z)\Cn. 

The validity of theorem^ can be formally justified by the following 
observation. For z close to A, Spec A{z) is the spectrum of A "bent" 
to the right half plane by the map w ^ {w — z)"^. Since the numerical 
ranges of A„(z) are far from zero (cf. lemma Ej), ||A„(z)~^|| are uni- 
formly bounded for all such z. When we add K{z), the truncation are 
not sectorial any longer, but since K{z) is of finite rank, for n large, 
An{z) + Kn{z) are still invertible, except for a few isolated points which 
correspond to perturbations of the isolated eigenvalue A of A{z) + K{z). 
When some of these isolated points is equal to zero, z G Spec2 (M|£„). 

Lemma 5. Let 6 > be small enough. Then there exist non-negative 
constants a, j3 and h, such that 

NumA„(2) C Conv [{[a , 00) + i[- ^ , I3]f] C {Re (2) > 6 > 0} 

for all n eN and z G {1^; — A| < S}. 

Proof. By substituting A for M in (0]), clearly l^umAn{z) C NnmA{z) 
for all n G N. Since A = A*, then A{z) is a normal operator for all 
2; G C, so that 

Num C Conv [Spec {A - z)^] z e C. 

Since A ^ Spec A, we can find 6 > small enough, such that — A| < 
6} does not intersect Spec A. Then there are a, (3 > 0, such that 

Spec (A-z) C ((-00, -a] + i[-p, /?]) U {[a, 00) + i[-p, /?]) 

for l^; — A| < 6. Hence 

Spec {A-z^ C {[a,oo) + i[-(3,(3]y. 

The b > can be found by choosing 6 > small enough, such that P 
is also small. 
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Lemma 6. Let ^ ip E DomH be such that Hip = vip and HnPn'ip 
uip. If H and Hn are invertible, and ||-f^^^|| < c where c > is 
independent of n, then 

as n oo. 
Proof 

\\H~^Pni^ - H-^i;\\ < c\\Pn1p - U-'HnPrM ^ 

as 72 — * OO. 

It is well known that for any bounded operator H and w ^ Num H, 
\\{H-w)-^\\ <dist {w,mmHy^ 

(cf. for instance [HI p.268]). Then according to lemma El for all n G N 

and l^; — A| < S, An{z) are invertible and there exists ci > independent 
of n and z, such that < ci. Since £ is an eigenspace of A{z), 

(H) ensures that An{z)Pnip A{z)ip as n ^ oo for all tp E S. Since 
z ^ Spec A, A{z) is invertible. Therefore lemma IHl yields 

(7) \\An{z)-'Pntp-A{z)-',p\\^0, ^eS 

for all \z — \\ < 6. 

The various constants Cj > that appear below are independent of 
z and n. 

In order to prove theorem^ we show that for 6 > small enough, we 
can find constants c > and N > uniform in z, such that an{z)~^ < c 
for all l-z — A| = S and n > N. In virtue of © and by applying lemma |21 
the above ensures the existence of > such that we can find 

Zn e {\z-X\ < 5}nSpec2 (M|£„) 

for all n > N. 

Let 5 > be small enough such that lemma El holds and 
SpecMn {|;z- A| < 6} = {A}. 

Denote by 

Ns:={zeC : \z - X\ =6}. 
Since A{z) + K{z) = (M - zY is invertible, for all z G N^, 

I + A{zy^K{z) = A{zy\A{z) + K{z)) 

is also invertible. Since the function ||/ + A{-y^K{-)\\ : Ns — > (0, oo) 
is continuous, it should be uniformly bounded. Put C2 > such that 

(8) \\{I + A{z)K{z))-^\\<c:^\ zeNs. 
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Let e > and z & Ns. Since Ns is compact, there exists a finite 
set of points {wj} C Ng satisfying the following property. Given any 
z E Ns, there is w G {wj} such that 

(9) \\A{zy^ - A{wy^\\ < £ and \w - z\ < e. 

Since dim£^ < oo and {wj} is finite, in virtue of there exists > 
such that 

\\A„,{w)-'P^^-A{wr'^\\<e, 
for all n > N, w E (wj) and ip E £. Thus, if w satisfies 

||A„(z)-ip„^-A(^)-Vll 

+ \\A{w)-'^-A{zr'n 

< \\An{zY^Pni^ - An{w)-^PrM + 2^ 

= \\A^{zy'A^{wy\A„{w) - A„(z))P„V^|| + 2£ 

< c\\\{w - z)P^[{w + z)- 2A]P^^\\ + 2e 
<\w — z\c^c^ + 2e < c^e. 

The existence of C3 is ensured by (H). By choosing e small enough, we 
can find > independent of z, such that 

\\An{z)-'PnK{z) - A{z)-'K{z)\\ < C2/2, n>N. 

Let X eH. In virtue of (jH)), 

C2\\PnX\\ < \\{I + A{z)-'K{z))PnX\\ 

< \\{I + Ar,{z)-'PnK{z))P^x\\ + 

+ \\Ar,{z)-'PnK{z) - A{z)-'K{z)\\\\PnX\\. 

Hence, for all n > N, 

(C2/2)||P„X|| < \\{I + An{z)-'PnK{z))P^x\\ 

< \\A^{Z)-'\\\\{A^{Z) + PnK{z)Pn)PnX\\ 

< Ci\\Pn{M - zYlCnPnXl 

Therefore Pn{M — z)'^\Cn is invertible and 

a^{z)-' = ||(P„(M-^)2|£„)-i|| < 2ci/c2 
when z E Ns and n> N. This ensures the validity of theorem ^ 
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4. An example of limiting set 



In this section we find the uniform hmiting set of second order spectra 
for a very simple model. We choose this example so that the spectral 
pollution is maximal when we apply the linear method. We also inves- 
tigate rank one perturbations of this model. These ideas are close to 
example I of [7| and those in [5]. 

Let E ^ (— TT, tt] be a finite union of semi-open intervals (0^,6^]. 
Denote by E'^ := {—71,11] \ E. Let Mf{x) := m{x)f{x) be the linear 
operator of multiplication by the symbol 



acting onH = L\-7t,7t). Let := Span {6"^'^^ . . . , e^"^'}. Then 
Mn is the (2n + 1) x (2n + 1) Toeplitz matrix {Mn)j,k = ^{j — k), 
where m denotes the Fourier coefficients of m. Notice that SpecM = 
SpeCessM = {±1}. 

Denote by T{m) the Toeplitz operator acting on whose symbol 
is m. By invoking Gohberg's theorem for piecewise continuous symbols 
(cf., for instance jU theorem 1.23]), we realize that Spec T(m) = [—1, 1]. 
The truncation of T(m) to the subspace 



where 6k{n) = 5k,n is the Kronecker delta symbol, equals M„ above. 
By virtue of a result due to Szego (cf. [H theorem 5.14]), for each 
/i G [—1, 1] there is a sequence /i„ G SpecM„ such that /i„ ^ A as n ^ 
00. Thus, if we choose the linear method with as approximating 
subspaces, each point in the spectral gap (—1, 1) of M is of spectral 
pollution. 

We now find the limiting set A(M). According to arguments in the 
proof of proposition 121 Spec2 (M|£„) lies in T for all n. We show that 
each C = e*^, — tt < 6' < tt, is in A. Indeed, by virtue of Gohberg's 
theorem for piecewise continuous symbols, the spectrum of T(m — C)^, 
the Toeplitz operator associated to the symbol (m — Q"^, equals the 
segment 




Cn = Span {6 



(N) 




P„(M - z^lCn = Q„T(m - Cf\C 



Since the essential spectrum of T{m — C,Y is a set with no interior, then 
each point in it is a Weyl's approximated eigenvalue (cf. theorem 
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10.10]). In particular so is the origin. Hence for all A; > 0, there is 
ijjk e P{N), \\ipk\\ = 1, such that \\T{m - OVfcll < 1/^- Since 

C < 11^ , II ' ^ "^-C m\ 

as n oo. Then C G S. Consequently proposition El ensures that 
A(M) = T. 

Let us discuss now how the points of Speca (M|£„) distribute along 
T. We saw in section 2 that 

z e Spec2 (M|£„) iff — ^ Spec M„. 

By construction, M„ coincides with the truncation of T(m) to the 
subspace Let 

{4"^ }Li := Spec2 {M\Cn) n {Im^ > 0} 

and 

(n) _(n) 

^("):=^^^G(-1,1). 
Then the mean of the points in Specg (M|£„) is 

n n 
-L \^ (n) . _(n) -L \^ (n) ^^rp 

— > 4 + 4 = - / ^A- = -Tr Mn 

k=l k=l 

= miO) = ^i\E\ - \E^\) 

for all n. Here | ■ | denotes the Lebesgue measure. 

A version of Szego's first limit theorem for Toeplitz operators 
theorem 5.10], permit us to say even more about the localization of 
Specj {M\Cn) in the limit n oo. According to [H corollary 5.12], for 
each Borel set B C [—1,1], 

1 ^ ^ \{x : m{x) E B}\ 

Then for all < £ < 1, 

10 - > 1 ^ - — - and - > 1 ^ ^ 

^ ' n 27r n ^ 27r 

«,("'6[-l,-l+e] 4"'e[l-e,l] 

as n — ^ oo. Consequently, although each point in T is in the limiting set 
A, for n large almost 100-|£^|/(27r) percent of the points in Speca (M|£„) 
cluster near 1 and the other 100 ■ |£"^|/(2'7r) percent cluster near —1. 
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A 


n 


lie Zn — 


Imzn 




ReZn + 


Imzn 




Re Zn — X 


-0.61803398 


85 


-0.64711164 


-0.59156988 


0.00130677 




120 


-0.64232258 


-0.59559824 


0.00092642 




155 


-0.63930167 


-0.59820046 


0.00071708 




190 


-0.63717720 


-0.60006016 


0.00058469 




225 


-0.63557976 


-0.60147516 


0.00049347 


1.61803398 


85 


1.58929960 


1.6448 


1953 


0.00097441 




120 


1.5939 


8716 


1.64069975 


0.00069052 




155 


1.59695260 


1.63804631 


0.00053453 




190 


1.59904216 


1.63615390 


0.00043595 




225 


1.60061565 


1.63471625 


0.00036803 



Table 1. Estimation of the isolated eigenvalues of M + 
K for a = 1, tp{x) = 1 and E = (0, vr]. The first column 
corresponds to the theoretical value. 



In order to envisage the numerical scope of theorem ^ we can test it 
against numerical data. For this we consider rank one perturbations of 
the above model. Let be a fixed vector such that \\ip\\ = 1, let a > 
and let Kf := a{f,ip)ip for all f ETC. We first compute explicitly the 
discrete spectrum oi M + K. 

Lemma 7. Let H he self-adjoint and such that Spec H = SpeCessH. 
Then the isolated eigenvalue of finite multiplicity of H + K are non- 
degenerate and they are the non-zero solutions A of 

(11) {{X-Hr'i,,^) = a-\ 

Proof U{H + K-X)(f) = Oioi (f)^0 and A ^ Spec H, then 

{H - \)(p + a{(p,ip)'ip = 0. 

If we normalize by {(f), ip) = a^^, then = (A — H)~'^ip. A substitution 
in the normalizing identity yields (fTT|) . Conversely suppose that (fTT|) 
holds. By putting := a(A — H)'^(f), we achieve {H + K — A)0 = 0. 

Hence, A 7^ ±1 is an eigenvalue of M + K iff 

^ ^ A-1 + A + 1 ""^ ' 

where = {2ti)-^ \^\^ dx. 

In tables 1 and 2 we show some numerical outputs produced by using 
the method of second order spectra applied to M -|- . We found Zn by 
adopting the algorithms for computing second order spectra available 
at the internet address mentioned in (7j. This value corresponds to 
the point in Specg (M + K\Cn), closest to the isolated non-degenerate 
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A 


n 


Re^;^ — 


Imzn 




ReZn + 


Imzn 




Re Zn — X 


-0.97901994 


85 


-0.99169545 


-0.9738 


4630 


0.00375093 




120 


-0.98897219 


-0.97406728 


0.00249979 




155 


-0.98740174 


-0.97435952 


0.00186068 




190 


-0.98635662 


-0.97465104 


0.00148388 




225 


-0.98561863 


-0.97491625 


0.00124750 


1.97901994 


85 


1.95326913 


2.00377483 


0.00049795 




120 


1.95700470 


2.00030477 


0.00036520 




155 


1.95961155 


1.99784090 


0.00029371 




190 


1.96164380 


1.99591836 


0.00023886 




225 


1.96314714 


1.99449943 


0.00019665 



Table 2. Estimation of the isolated eigenvalues of M + 
K for a = 1, iplx) = 1 and E = (— 157r/16, vr]. The first 
column corresponds to the theoretical value. 



eigenvalue A. We found the theoretical value of A by using (fT^ . In 
both tables it is remarkable that although bound (j21) only provides 
estimation for the first digit of A at = 225, the actual value of Re Zn 
is correct up to three digits. For table 2 we chose an extreme case 
where one of the non-degenerate eigenvalues is close to the point -1 of 
the essential spectrum. It is remarkable that the even the step n = 85, 
which only takes a few second to run in a PC, already detects the 
presence of this eigenvalue. 

Figures 1, 2 and 3, shows how the rest of the second order spectrum 
in these two cases distribute along the complex plane for three different 
values of n. Notice that the clustering predicted for the unperturbed 
case seems to be largely kept with the exception of some few points 
that approximate the isolated eigenvalues. 

5. The condition (H) 

As in the linear case, an effective numerical implementation of the 
method of second order spectra to unbounded operators, demands veri- 
fying (H) without having much information about the eigenf unctions ip. 
This is achieved sometimes by means of a perturbative argumentation. 
We do not claim that this is easy in general, but in some situations the 
following standard argument can be useful. 

The general strategy involves a dominating operator X, which has 
compact resolvent and a complete set of eigenfunctions that can be 
found explicitly. Both operators M'^, k = 1,2, should be relatively 
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X-bounded, in the sense that DomX C DomM^ and 

IIM'^mII < afellXwII + bk\\ u\\, M G DomX, 

for uniform constants > and bk > 0. The best such that the 
above holds for some bk, is called relative bound. Put as approximating 
sequence {4>n}'^=i, the set of eigenf unctions of X. If the eigenfunction 
ip of M associated to the eigenvalue A lies in DomX, by virtue of the 
spectral theorem, X and P„ commute, and hence 

||xp„^ - x^ii = \\PnX^ - x^pw ^ 

as n oo. Then WM'^Pnip — M''ip\\ 0, ensuring (H). 

The following result illustrates this strategy. Recall that the class 
Ki consists of all : M — > M such that 

sup / \V{y)\dy < oo. 

xm J\x-y\<l 

Theorem 8. Let 0„ be the eigenf unctions of the harmonic oscillator 
{-dl + x^) acting on L'^{M). LetV -.M. — ^ M 6e such that 

V e[Kir} W^2'^(M)] + iy°°'2(M). 

Then (H) holds for every eigenvalue and eigenfunction of the operator 
M = -dl + V acting on ^^(M). 

Analogous results can be found for higher dimensions. 

Proof. Since V G L'^ + L°°, multiplication by V is relatively 9^-bounded 
with bound (cf. |3j §1.2]), hence M is relatively 9^-bounded with 
bound 0. 

We prove that is relatively 9^-bounded with bound 1. For this, 
notice that 

{-dl + V)\ = dlu - 2Vdlu - 2V'd^u + (V^ - V")u. 

Since V E L"^ + L°°, multiplication by V is (9^-bounded with relative 
bound and hence 2Vdl is 9^-bounded with relative bound 0. Since 
G + L°°, a similar reasoning ensures that 2V'dx is 9^-bounded 
with relative bound 0. Since V" lies in + L°°, multiplication by V" 
is 9^-bounded with relative bound 0. Since V'^ is 9^1^-bounded with 
relative bound and 

dlVu = V"u + 2V'dxU + Vdlu, 

then V"^ is 9^-bounded with relative bound 0. This ensures that is 
9^-bounded with relative bound 1. 
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By virtue of Kato-Rellich theorem, above we do not have to worry 
about the domains. Here we are considering 

DomM = 1^2,2 DomM2 = 1^2,4 

The mentioned resuh ensures that M (and hence M^) are self-adjoint 
in these domains. 

Put X = {—d1 + x^y. Since V G Ki, then \ii{x)\ < ce"'"'''! for some 
constants a > and c > (cf. [HJ §C.3]). Together with the inclusion 
V' e IV^'^, this bound ensures ip G Dom (—9^ + x^) = W^^'^ fl Dom (x^) 
and (—'?/'" + x^?/') G VT^'^ fl Dom(x^). The only non-trivial facts of 
the latter assertion are, perhaps, the inclusions x'^ip G W'^''^ and ip" G 
Dom(x^). The first follows from the second, by differentiating twice 
the term x'^ip and noticing that {xip')' G L^. The second inclusion is 
achieved by means of the following trick. Since Mip is eigenvector of 
M, I - <(x) + V{x)ij{x)\ < ce-"l^l. Then 

llxVll < \\x^{-ip" + 1/^)11 + llxV^II 

<ci+ x^|V(x)nV^(x)|Mx^ ^ 

1/2 



The latter integral is bounded because of ^ G + L°°, therefore 
tp" G Dom(x^). Thus G DomX. 

Finally we show that is (9^ + x^) ^-bounded. For this, let 

A := 2-^/\x + d^) and A* = 2-1/^(3, _ 

Then (-5^ + x^) = 2(A/l* - 1) and 9^ = (A - Thus the desired 

property follows from the identity (cf. jlll eq.(X.28)]) 

\\A*'---A*"u\\ < c\\{-dl + x^y'/^l n = l,2,... 

where A**" is either A or A*. This is easily shown by induction and 
using the estimate 

U-dl + x^)''/\\\ < ||(-9^ + x2)"/2m||, k<n. 

This completes the proof of the theorem. 

In the above result, we have in mind the case V = W + S, where W 
is periodic and S is in L^(]R) + L°°(]R)e. Under this hypothesis, it is 
well known that the essential spectrum of M has a band gap structure 
determined solely by W, whereas the perturbation S can produce non- 
empty discrete spectrum (cf. e.g. [121 Chapter XIII. 16] and also ;IU]). 
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Figure 1. Spec2 {M + K\Cn) for three values of n. The 
operator corresponds to a = 1, ip{x) = 1 and E = (0, tt]. 
See table 1. 
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Figure 2. Spec2 (M + K\Cn) for three values of n. 
The operator corresponds to a = 1, ip{x) = 1 and 
E = (-157r/16,7r]. See table 2. 
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Figure 3. Two zoom pictures of figure 2. Clustering 
near ±1. The star shows the location of the theoretical 
eigenvalue. 



